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ABSTRACT 

We introduce the numbers of hot and cold spots, n h and n c , of excursion 
sets of the CMB temperature anisotropy maps as statistical observables that can 
discriminate different non-Gaussian models. We numerically compute them from 
simulations of non-Gaussian CMB temperature fluctuation maps. The first kind 
of non-Gaussian model we study is the local type primordial non-Gaussianity. 
The second kind of models have some specific form of the probability distribution 
function from which the temperature fluctuation value at each pixel is drawn, 
obtained using HEALPIX. We find the characteristic non-Gaussian deviation 
shapes of rih and n c , which is distinct for each of the models under consideration. 
We further demonstrate that nh and n c carry additional information compared to 
the genus, which is just their linear combination, making them valuable additions 
to the Minkowski Functionals in constraining non-Gaussianity. 

Subject headings: CMB, non-Gaussianity 



1. Introduction 

The statistical nature of the cosmic microwave background (CMB) radiation temper- 
ature fluctuation field that we see today must be predominantly inherited from those of 
the primordial density fluctuations. Inflation is currently the forerunner amongst possible 
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mechanisms that could have produced the primordial density perturbations. All models of 
inflation, in general, predict some amount of deviation of these fluctuations from a Gaus- 
sian distribution, with the details of the deviations being model dependent. The knowledge 
of these deviations, if observed, for example in the CMB, will thus be of much value in 
distinguishing between various models of inflation. The observational search for primor- 
dial non-Gaussianity, however, is not easy since various observational effects can mask the 
true CMB signal. Given this difficulty the need for efficient, sensitive and complementary 
observables that can characterize non-Gaussian deviations cannot be over emphasized. 

Popular statistical measures of non-Gaussianity that can be obtained from the CMB 
fall un der two broad categories. The first are observables that are defined in harmonic 



space ( IKomatsu et al.ll201ll ; ISmidt et al.ll2010t ICurto et al.ll201ll ; IVielva fc Sanzll2010l ) such 
as the bispectrum,trispectrum, wavelets, the spherical Mexican hat wavelet etc. The sec- 
ond category consists of those that directly exploit the geometric and topological proper- 
ties of the tempera ture fluctuation field. A popular class of observables, the Minkowski 



Functionals (MFs) (ITomita Il986t IColeslll988 



Gott et el1.M199CH : ISchmalzing & Gorskilll998 



Winitzki fc Kosowskyl Il998[ ). have long been applied to constrain non-Gaussianity in the 
CMB. Considerable progress has been made i n understanding them analytically for weakly 



non- Gaussian random fields ( Matsubara 



2003 



Pogosvan et all 120091 ; lMatsubarall201d : iGav et all 120121 ). Of the three MFs that can be de- 



Hikage et al.ll2006l . 120081 : IKomatsu et al.ll2009 



fined for a 2-dimensional random field, the third one known as the genus is a topological 
quantity which depend on the global properties of the random field. It is given by the differ- 
ence of the numbers of hot and cold spots at any given temperature fluctuation field value. 
The genus and other MFs have non-Gaussian deviation shapes which is characteristic of the 
non-Gaussian model. The non-Gaussian deviation shape of each observable tells us what 
field values are best probed by the observable, and these are the values where the deviations 
are the largest. 

In this paper, we introduce the numbers of hot and cold spots as statistical observ- 
ables in their own right. Just like the genus they are topological quantities which depend 
only on the global properties of the temperature fluc tuation field. They have been studied 
earlier by Coles and Barrow (jColes fc Barrow! 119871 ) and approximate formula are known 
for Gaussian and some non-Gaussian random fields JVanmarcke!ll983h . Our goal here is 
to determine their sensitivity and non-Gaussian deviation shapes as signatures of various 
kinds of non-Gaussianities. We have computed them using numerical methods from simu- 
lated non-Gaussian CMB maps. We first compute them for simulatio ns containing the so 
called loc al type primordial non-Gaus sian model parametrized by f^L (jLiguori et al.l 120031 ) 
and Qnl ( jChingangbam fc Park! 120091) . We also compute them using non- Gaussian simu- 
lations provided by HEALPIX (IRocha et ajj 120051 ; iGorski et al.l 120051 ). We have obtained 
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the characteristic non-Gaussian deviations for these different types of input primordial non- 
Gaussianity. 

We do not intend to address all the issues of observational systematic effects in this first 
paper and focus on the theoretical understanding of their behaviour upon potential non- 
Gaussianity. The expectation in introducing them is that we can get additional information 
about non-Gaussian fields by using these two observables in addition to the genus and other 
MFs. Since the genus is given by the difference between these two numbers, in taking 
their linear combination we are throwing away some information. This expectation will be 
most justified if the number of hot spots is independent of the number of cold spots. As 
we will see in section 3, this is not always the case and whether they are independent or 
not depends on the non-Gaussian model. For example, for the local type primordial non- 
Gaussian model parametrized by /jvl and g^L, we find that they are related to each other in 
a specific way. However, even for the models where they are related, there is still additional 
gain of information coming from the fact that they have non-Gaussian deviation shapes 
which are quite distinct from those of the genus and the other MFs and hence they best 
extract non-Gaussian deviations of the field at values different from the MFs. Therefore, 
they carry inormation complementary to the MFs. We further calculate the uncertainties 
in the numbers of hot and cold spots, taking into account cosmic variance, instrumental 
noise and sample boundaries at our choice of smoothing scales. We demonstrate that there 
exists additional information in the numbers of hot and cold spots compared to the genus, 
as shown in section 4. 

This paper is organized as follows: in section 2 we briefly describe excursion sets and hot 
and cold spots, followed by our method for numerically computing them. We then show the 
results for the numbers of hot and cold spots computed from Gaussian CMB temperature 
fluctuation maps. In section 3 we present the non-Gaussian deviations of the numbers of 
hot and cold spots for the different kinds of non-Gaussianity that we have studied, first, the 
local type non-Gaussianity and secondly for the input non-Gaussian PDF models obtained 
using HEALPIX. In section 4 we analyze how observational effects such as beam profiles, 
instrument noise and incomplete sky coverage, affect the numbers of hot and cold spots 
by computing them from simulations to which these effects have been added. Further, we 
compare their sensitivity to non-Gaussianity with that of the genus and show that they 
can give more information than the genus. We end with a summary and discussion of the 
implications of our results in section 5. 
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2. Hot and cold spots counts 



Let / = (T(n)—T )/T denote the CMB temperature anisotropy field, where, T denotes 
the mean temperature. By rescaling / by its rms value, <To, we can define the threshold 
temperature, v = f/a . At each value of u, if we consider the set of all pixels that have 
values equal to or above v we obtain what is usually referred to as an excursion set. This set 
consists of many connected regions into which the temperature field 'manifold' has fractured, 
and holes within those regions due to the excluded pixels. As is commonly done in the 
literature, we call each connected region a hot spot and each hole a cold spot. For an excursion 
set indexed by u, we define: 

• rih = number of hot spots, and 

• n c = number of cold spots. 

As we change v, the excursion sets behave as though they are a one-parameter family of 
spaces parametrized by u, and their properties such as the numbers of the hot and cold 
spots, change systematically as a function of v. 

We can relate the numbers of hot and cold spots to the numbers of closed iso-temperature 
contours. The boundaries of each excursion set are iso-temperature contours for the corre- 
sponding v. We can assign an orientation to each of the contours - positive for those that 
enclose hot spots and negative for the ones that enclose cold spots, rih and n c are then simply 
counts of closed positive and negative orientation contours, respectively. For the purpose 
of illustration, in Fig. ([1]), we have shown a smoothed Gaussian random field defined on a 
square with periodic boundary condition. The left panel shows the full field. The middle 
panel shows the excursion set for the same field for v = 0. Connected regions and holes 
are clearly visible. The right panel shows the boundary or iso-temperature contours for the 
same excursion set, red enclosing hot spots and blue enclosing cold spots. 

Mathematically, we can express rih and n c as line integrals 



where K is the total curvature of iso-temperature contours for each v. C + denotes contours 
that enclose hot spots while C_ denotes contours that enclose cold spots. The genus, g, is 
given by a linear combination of and n c : 






n h {v) - n c {y). 



(2) 
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Fig. 1. — Left panel: A patch of a smoothed Gaussian fluctuation field. The field is defined 
on a square with periodic boundary condition. So field values that are at opposite ends of 
the square are identified. The y-axis gives the level value, v. Middle panel: The red regions 
gives the excursion set for v = 0. The set is fragmented into several contiguous or connected 
regions. Some of the connected regions have holes within them. Each connected region is 
called a hot spot while each hole is called a cold spot. Right panel: Iso-temperature contours 
enclosing the excursion region (red lines) and holes (blue lines) for v — 0. One can see 
partial contours that are located at extreme ends of the square which together form closed 
contours, is the number of isolated connected regions, which can be obtained by counting 
the closed red contours. n c is the number of holes within the connected regions and can be 
obtained by counting the blue contours. 

For a Gaussian fluctuation field, g is given by the expression 

g{v) = Ave-^l\ (3) 

where the amplitude A carries the physical information about the field and is given by 
A = (1/27T 3 ) (ai/ao) 2 , with cxi being the rms of the gradient field. 

Note that we could equivalently define our excursion set at each v to consist of the pixels 
having values below v, in which case the definitions of and n c would get interchanged. 



2.1. Numerical method for computing the numbers of hot and cold spots 



Our method for computing the numbers of hot and cold spots is based on the method 
for computing the genus outlined in (IGott et al.lll990l ). It is essentially an implementation of 
Eq. ([T]) based on connecting iso-temperature pixels with the information of the orientation 
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of the contour retained. At the end contours with the same orientation are counted to get 
rih and n c . 
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Fig. 2. — rih, n c and g measured from Gaussian simulations smoothed with FWHM=30'. 
The y axis values are per unit area of the sphere. The plots are average over 200 simulations. 



In Fig. (j2j) we have shown rih, n c and g versus v obtained by averaging over measure- 
ments from 200 simulated Gaussian CMB anisotropy maps with the HEALPIX resolution 
parameter N 'side = 512 and smoothed with a Gaussian filter with FWHM= 30'. The simu- 



lations have ACDM parameter values given by WMAP 5 years data (IKomatsu et al.ll2009l ). 
The Gaussian genus formula serves to test the accuracy of the numerical computation of rih 
and n c (unless of course there is some error whic h contributes equ ally to both rih and n c and 
cancels out for the genus). It has been shown in (IPark et al.ll2012l ) that precise details of the 
numbers of hot and cold spots in terms of Betti numbers for Gaussian random fields such as 
the amplitude and location of peaks vary significantly as we vary the power index, n, of the 
input three-dimensional power spectrum P(k) oc k~ n . The trend is that as n increases the 
amplitude increases and the peak shifts closer to v — 0. For n ~ 3, which is relevant for the 
CMB, the result that we have obtained is in agreement with this trend. 



3. Non-Gaussian deviations of the numbers of hot and cold spots 



In this section we compute the non- Gaussian deviations of numbers of hot and cold 
spots for different models of primordial non-Gaussianity. 
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3.1. Local type primordial non-Gaussianity 



We consider the following expansion to cubic order of the primordial gravitational po- 
tential: 



$(*) = + f NL (($ G (f)) 2 - (($ G ) 2 » + g NL ($ G (x)) 3 + . . . , (4) 

where $ G is a Gaussian potential and /jvl and g^L are constants which parametrize the first 
and second order non-linearities respectively, in the gravitational potential. Then, expanding 
the CMB temperature fluctuation field in multipoles, as f — J2 em ae m Yi m , we can calculate 
o-im by convolving $ with the full radiation transfer function A^, as 

d 3 k 



4tt( 



§{k,ti) A £ (k,t ) Yf m (k) 



(5) 



We u se simulations of non-Gaussian CMB maps (ILiguori et al.l 12003c IChingangbam fc Park 
20091 ) which have Eq. (j4j) as the input potential to calculate ai m . The input power spectrum 
of $ G is given as P$(k) = (A n /k 3 )(k/ko) ri3 ~ 1 , where A , n s and k are taken from WMAP 



5 years ACDM parameter values (IKomatsu et al.l 120091 ) . The simulation resolution used is 
given by NSIDE = 512, as in section (12. ip. We use Ai calculated with all perturbation 
terms kept to linear order (ISeljak fc Zaldarriagalll996l ) and hence the non-Gaussianity that 
shows up in the resulting CMB temperature field is a direct consequence of the primordial 
input. We have computed rih and n c for three kinds of simulations - pure /jyx, pure g^L and 
a mixture of the two. In order to quantify the non-Gaussian deviations we define 



nf G 



n, 



G 



(6) 



where % stands for h or c, the index G stands for Gaussian and NG for non-Gaussian. Plots 
are shown normalized by n G,max , which is the amplitude of nf. 

In Fig. ([3]) we have plotted An^ and An c versus v for pure f^L and pure g^L cases. We 
have used the values /jvx = ±100 and g^i = ±1 x 10 6 and shown plots for two smoothing 
angles - FWHM= 30' and 90'. For each case, it is important to note that Arii has a 
characteristic non-Gaussian deviation shape and they can be easily distinguished from each 
other. There is slight variation of the deviation shapes as functions of the smoothing angle. 
Roughly speaking, the magnitude of the deviation at higher threshold values \u\ > 2 are 
larger for larger smoothing angles. An interesting observation is that for each case we can 
see that and n c are correlated as 



An h (vJ NL ) 
An h (v, g NL ) 



-An c 
AnJ- 



'-v, Inl) 



(7) 



In Fig. (TjJ we have plotted Anh versus An c at each value of v for pure /nl (green, solid 
line) and for pure g^L (brown, dotted) models. This is another way of visualizing the 
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characteristics of the non-Gaussianity caused by f^L or g^i- Fig. © shows Arii for the case 
when both f NL and g^L contribute to the primordial non-Gaussianity for the same smoothing 
angles and parameter values. If we compare with non-Gaussia n deviations of the genus (see 
Fig. (2) of ( IHikage et al.l 120061 ) for pure f NL case and Fig. (4) of (jChingangbam k, Parkl|2009l ) 
for pure g^L case), we find that the amplitude of the deviations of nh and n c are smaller by 
about factor of two. 

We can get some idea about the dependence of Arih and An c on f NT , and Qnl from 



the analytic exp ressions of the non-Gaussian deviation of the genus (IHikage et al.l 12006 



Matsubaral 120031 ) . which is given as an expansion in powers of cr . For pure /jvl and pure 
gNL cases, keeping the genus expansion upto cr and (Tq orders, respectively, the genus non- 
Gaussian deviations have linear dependence on f NL and g^L- When both f NL and g^L are 
present , then at <Tq order there must be cross terms containing both f NL and g^L- Hence 
the non-Gaussian deviation of the genus will not be a simple linear combination of deviation 
terms depending on f^L and g^L independently. Since the genus is just the subtraction of 
n c from nh, we can expect n c and to behave in a roughly similar fashion. 



3.2. HEALPIX non-Gaussian models 



We have generated n on-Gaussian maps using the HEALPIX routine sky_ng_sim ( Mocha et al 



20051 ; iGorski et al.ll2005l ). This program implements two kinds of non-Gaussian models. The 
first is a model where the input probability distribution function is taken to be an expansion 
in excitepd states of the simple harmonic oscillator (SHO model), as given below: 

2 



P{f) 



( 



./ 



where Hi are Hermite polynomials, Cj are normalization constants, <Jq is the variance of 
the Gaussian PDF and «j, for % > 1, are free parameters. «o is constrained to be «o = 
yT— Yl~i \ a i\ 2 - F° r our simulations we have kept terms upto n = 2 such that ot\ and «2 are 
non-zero. The second non-Gaussian model has the input PDF of the temperature field as an 
even power of a Gaussian PDF (Gaussian power model), with the temperature fluctuation 
value of the k— th pixel given by 

f(k) = g 2p (k), 

where g is a zero mean, unit variance Gaussian variable, and p is chosen to be a positive 
integer. We have used p = 1 for our simulations. 

Fig. ([6]) shows the non-Gaussian deviations of and n c for these two models. Arij is 
again defined as given in Eq. ([6]). The left panel shows the deviations for SHO model and 
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right panel shows for Gaussian power model. As in the local non-Gaussianity case we can 
see characteristic deviations for each type of non-Gaussianity. 

For these models they do not have any simple correlation between An^ and An c such as 
what we have seen for the local type non-Gaussianity indicating that they carry information 
independent from each other. Fig. (JTj) shows An^ versus An c at each value of v for these 
two models. These plots are a different way to characterize the type of non-Gaussianity of 
the HEALPIX non-Gaussian models. 

4. Statistical sensitivity of the numbers of hot and cold spots to f^i and g^i 

To analyze the statistical power of nn and n c in realistic situations we measure them 
from simulations to which observational effects have been added. The observational effects 
are pixel window function, beam profile for each differential assembly (DA) and Gaussian 
noise realizations for each DA that follow the noise pattern, followed by Galaxy and point 
source masking. We then coadd Q, V and W DA's with appropriate weights obtained fron 
the inverse of the full-sky averaged pixel-noise variance in each DA, and then smooth the 
field. For the Galaxy masking we use the KQ75 mask. In Fig. (jSJ) we have shown the 
sample variance error bars obtained from the 200 Q + V + W coadded maps prepared as 
described above. It is immediately noticeable that the error bars for g^L is larger than those 
of /jvl at each corresponding smoothing angle. This can be understood from Eq. (J4]) as 
follows. Suppose we have a perfectly Gaussian field <p G and another 'Gaussian' field with 
slight statistical fluctuations cf>' = 4> G (1 + D), where D quantifies the fluctuation. Then the 
deviation from the Gaussian field when f NL and g^L contributions are present is given by 
A4> NG = D(p G + 2D f NL (4> G ) 2 + 3D g NL (<f) G ) 3 . Hence statistical fluctuations seen for g NL will 
be larger than those for fj^L- 

As a simple way of estimating the statistical disciminating power of the numbers of 
hot and cold spots in comparision to the genus we integrate the absolute values of the non- 
Gaussian deviations measured in units of the corresponding sample variances from v = —3 
to 3. Let us denote it by A. For M threshold levels with spacing Au, we can calculate it as, 

M 

A = AuJ2(\^0{i)\/0 G ^) /a s (i), (9) 

where O can be either g, or n c , and a s (i) are the respective sample variances at each 
threshold level i. For our case, M — 31 and Av = 0.2. The resulting values are shown 
in Table (JT|). We find considerably larger values of A for rih and n c compared to g both 
for /jvx and g^L type non-Gaussianities, at the smoothing angles we have considered. This 
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demonstrates that there is loss of statistical power for detecting the presence of non-Gaussian 
deviations when we combine and n c to get the genus. 



5. Conclusion 

We have introduced the numbers of hot and cold spots of the CMB temperature fluctua- 
tion field as statistical observables in their own right and propose to use them as discriminants 
of non-Gaussianity. We have studied the theoretical predictions for the numbers of hot and 
cold spots and their expected non-Gaussian deviations for various kinds of non-Gaussianities. 
We have calculated them using numerical methods from simulated CMB maps containing 
the different non-Gaussian models as inputs. The first type of input non-Gaussian model 
we studied is the so called local type primordial non-Gaussianity, parametrized by f^L and 
g^L at the first and second order non-linearity respectively, of the perturbative expansion 
of the primordial gravitational potential. This gravitational potential is convolved with the 
full radiation transfer kept to linear order to obtain the simulated CMB temparature field, 
and hence the non-Gaussian deviations seen in the numbers of hot and cold spots are direct 
probes of the primordial non-Gaussianity. For these local primordial non-Gaussian models, 
what we find is that rih and n c are correlated as given by Eq. (I7j). The strengths of the non- 
Gaussian deviations of nh, n c and g are large at different regions of v and hence each of them 
best probe regions of the field values specific to it. Therefore, they provide complementary 
information. Further, we have demonstrated that there exists additional information in the 
numbers of hot and cold spots compared to their linear combination given by the genus. 

The second class of non- Gaussian models that we have considered assumes specific forms 
of the PDF of the temperature fluctuation field. In particular, we studied a model where 
the simulated temperature value at each pixel is drawn from a PDF given as an expansion 
in simple harmonic oscillator states. We also studied another model where the temperature 
fluctuation values at each pixel is given as even powers of a number drawn from a Gaussian 
distributed field. Note that the assumption of the form of the PDF does not tell us anything 
about the physical source of the non-Gaussianity. Even though the physical origin of the 
non-Gaussianity is not clear, they are quite interesting models because they provide examples 
of non-Gaussian models where nh and n c are not correlated. 



It is interesting to compa re the shapes of the numbers of hot and cold spots with those of 
maxima and min ima counts (jPogosyan et al.ll201ll ; iBardeen et al.lll986t iBond &: Efstathiou 
19871 ; lAdlerl Il98ll ) . The shapes of and n c are roughly similar to the maxima and minima 
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counts, respectively, though the precise shape information such as peak location and the 
amplitude is quite different. Note that the extrema counts will measure more number of 
objects per unit area in comparision to the numbers of hot and cold spots since a typical 
connected/hole region can have more than one maxima/minima. As \u\ becomes much larger 
than one, rih and n c should tend towards the maxima and minima counts, respectively. 

Our next goal is to apply the number of hot and cold spots to observational data and 
constrain /tvl and g^L- If would also be very useful to have their analytic expressions. We 
are presently working towards these directions. 

We thank Korea Institute for Advanced Study for providing computing resources (KIAS 
Center for Advanced Computation Linux Cluster System QUEST) where the local non- 
Gaussian simulations used in this paper were computed. We also acknowledge use of the 
Hydra cluster at the Indian Institute of Astrophysics for a part of the analysis. We acknowl- 
edge use of the HEALPIX package. 

REFERENCES 

Adler, R. J. 1981, The Geometry of Random Fields, ed. Adler, R. J. 

Bardeen, J. M., Bond, J. R., Kaiser, N., & Szalay, A. S. 1986, ApJ, 304, 15 

Bond, J. R., & Efstathiou, G. 1987, MNRAS, 226, 655 

Chingangbam, P., & Park, C. 2009, J. Cosmology Astropart. Phys., 12, 19 

Coles, P. 1988, MNRAS, 234, 509 

Coles, P., & Barrow, J. D. 1987, MNRAS, 228, 407 

Curto, A., Martinez-Gonzalez, E., Barreiro, R. B., & Hobson, M. P. 2011, MNRAS, 417, 488 

Gay, C, Pichon, C, & Pogosyan, D. 2012, Phys. Rev. D, 85, 023011 

Gorski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759 

Gott, III, J. R., Park, C, Juszkiewicz, R., et al. 1990, ApJ, 352, 1 

Hikage, C, Komatsu, E., & Matsubara, T. 2006, ApJ, 653, 11 

Hikage, C, Matsubara, T., Coles, P., et al. 2008, MNRAS, 389, 1439 



- 12 - 



Komatsu, E., Dunkley, J., Nolta, M. R., et al. 2009, ApJS, 180, 330 
Komatsu, E., Smith, K. M., Dunkley, J., et al. 2011, ApJS, 192, 18 
Liguori, M., Matarrese, S., & Moscardini, L. 2003, ApJ, 597, 57 
Matsubara, T. 2003, ApJ, 584, 1 
-. 2010, Phys. Rev. D, 81, 083505 

Park, C, Pratyush, P., Chingangbam, P., et al. 2012, ApJ to be submitted 

Pogosyan, D., Gay, C, & Pichon, C. 2009, Phys. Rev. D, 80, 081301 

Pogosyan, D., Pichon, C, & Gay, C. 2011, Phys. Rev. D, 84, 083510 

Rocha, G., Hobson, M. P., Smith, S., Ferreira, P., & Challinor, A. 2005, MNRAS, 357, 1 

Schmalzing, J., & Gorski, K. M. 1998, MNRAS, 297, 355 

Seljak, U., & Zaldarriaga, M. 1996, ApJ, 469, 437 

Smidt, J., Amblard, A., Byrnes, C. T., et al. 2010, Phys. Rev. D, 81, 123007 
Tomita, H. 1986, Progress of Theoretical Physics, 76, 952 

Vanmarcke, E. 1983, Random fields: Analysis and synthesis, ed. Vanmarcke, E. 
Vielva, P., & Sanz, J. L. 2010, MNRAS, 404, 895 
Winitzki, S., k Kosowsky, A. 1998, New A, 3, 75 



This preprint was prepared with the A AS IATgX macros v5.2. 



13 



| 0.005 



a j3 



o 



< -0.005 



g 0.005 

R 



< -0.005 



FWH 

--' \, n 


M = 30' 
/ \/ 


1 








FWH 

^~^\^ ^ 


M = 90' 

/ V S 

/ \ / 




1 i'NL 

III 


\-' \ 

oo, g NL =0 
100 g NL = 

1 1 1 





-1 1 

v 



FWHM = 30' 




FWHM=90' 




-0.005 



f NL =100, g NL =0 

- f NL =iioo g NL =0 



-1 1 

V 



a 



6 





I 


FWHM = 30' 






0.005 






X 
cO 

s 


0.005 









o' O 





























0.005 






< 


-0.005 






FWHM = 90' 






0.005 






CO 


0.005 








E 











O o 

c 





















0.005 




f N L=0, g NL =10 6 


<a 


-0.005 




l 


_ f NL = g NL =-10° 

I I I I I I I I I I 









1/ 



III 

FWH 
~ N /~~ 


M = 30' 

^ ^ 

✓ 


FWH 
\ 


Vl=90' 


_ f NL =0 

I 1 I 1 I 


/ ^ _ — 

/ gNL=10 6 

' gNL T~ 10 t , - 





v 



Fig. 3. — Non-Gaussian deviations of rih and n c for pure f^L (upper panels) and pure gNL 
(lower panels) input primordial non-Gaussianity at two smoothing angles - FWHM= 30' 
and 90'. An^ is defined as given in Eq. fl6]). nf ,max is the maximum value of nf(v). The 
y axis values are per unit area of the sphere. The simulations have WMAP 5-years ACDM 
parameter values. The results are average over 200 simulations. 
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Fig. 4. — Plots of Ant versus An c for pure f^L (green, solid line) and for pure g^L (brown, 
dotted) models. This is a different way of showing the characteristics of the non-Gaussianity 
caused by the /nl and qnl terms. 
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Table 1: Table showing values of A defined in Eq. (jHJ) for g, nh and n c . 
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Fig. 5. — Same as in Fig. (EJ) for mixture of /nl and g^L input primordial non-Gaussianity. 
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Fig. 6. — Left panel: Arii for the SHO model. We have used a\ = 0.6 and cti = 0.6. Right 
panel: Arii for Gaussian power model. We have used p — 1. The plots are average over 50 
maps. 




Fig. 7. — Plots of Anh versus An c for SHO (green, solid line) and Gaussian power (brown, 
dotted line) models. 
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Fig. 8. — Non-Gaussian deviations of n.^ and n c for f^L (upper panels) and gNL (lower 
panels) calculated after adding observational effects, namely pixel window function, beam 
profiles, noise for each DA and galaxy and point source masking, to the simulations. These 
calculations are from Q + V + W coadded maps. The error bars are the sample variance 
from 200 simulations. 



